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Abstract 

In the two-dimensional framework, the surface gravity of a (classical) black hole is independent 
of its mass M. As a consequence, the Hawking temperature and outflux are also independent of M 
at the large-M limit. (This contrasts with the four-dimensional framework, in which the surface 
gravity and temperature scale as 1/M.) However, when the semiclassical backreaction effects on 
the black- hole geometry are taken into account, the surface gravity is no longer M-independent, 
and the same applies to the Hawking temperature and outflux. This effect, which vanishes at the 
large-M limit, increases with decreasing M. Here we analyze the semiclassical field equations for 
a two-dimensional static black hole, and calculate the leading-order backreaction effect (oc 1/M) 
on the Hawking temperature and outflux. We then confirm our analytical result by numerically 
integrating the semiclassical field equations. 



I. INTRODUCTION 



In classical General Relativity, a black hole (BH) is absolutely black: It does not emit any 
radiation. The situation changes, however, when quantum effects are taken into account. 
The semiclassical extension of General Relativity considers quantum fields which live on the 
background of a well-defined classical geometry (e.g. a black-hole). Within this framework 
it was found |lj that a black hole actually has a finite temperature, the Hawking temperature 
Th. Accordingly the BH emits thermal radiation, and evaporates within a finite time. 

Hawking's analysis pQ revealed that the temperature of the semiclassical BH is Th = 
k/2"7t, where k is the surface gravity of the BH. Throughout this paper we use General- 
Relativistic units c = G = 1 (and the same for the Boltzmann constant), and we also 
set h — 1 (following Ref. [2]). The temperature Th is thus uniquely determined by the 
background BH geometry. For a 4-dimensional (4D) Schwarzschild BH of mass M the 
surface gravity is k = 1/(4M), hence Tjj = l/(87rM). 

In the semiclassical theory, the quantum field yields an Energy-momentum contribution 
T a p, to be inserted at the right-hand side of the semiclassical Einstein equation G a p = 8nT a/ 3. 
This renormalized stress-energy tensor T a p is a tensor field in spacetime, which depends 
on the spacetime geometry (as well as on the quantum state of the matter field under 
consideration). The Hawking radiation is the most obvious manifestation of this semiclassical 
T a/3 : It is the outging component of T a/3 , evaluated at future null infinity (FNI). But obviously 
T a/3 includes other components as well, and is also position-dependent. For example, an 
evaporating BH must be endowed with a negative ingoing component of T a p(x^) at the 
horizon: It is this ingoing component which is directly responsible for the steady decrease 
of the BH mass. 

The classical Schwarzschild solution is a vacuum solution of the Einstein equation. Obvi- 
ously, the semiclassical contribution T a p must modify the BH geometry, which will no longer 
be a pure vacuum solution. Instead, the BH metric g a p is to be determined now from the 
semiclassical Einstein equation G a/ 3 = 8nT a /3. 

Since in a semiclassical BH g a p is no longer the Schwarzschild geometry, the BH's surface 
gravity will deviate from its classical value k = 1/(4M). This deviation is small for a 
macroscpic BH (namely, M 3> m p , where m p denotes the Planck mass), and is expected 
to decrease with increasing M. This change in k yields a corresponding change in the 
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BH temperature and outflux. The main objective of this paper is to explore this (mass- 
dependent) change in the temperature and outflux of a semiclassical BH, caused by the 
deviation of the background geometry from the classical one. (Though, we shall actually 
tackle this problem in two rather than four spacetime dimensions, for reasons which we 
shortly explian.) 

Presently there is no known explicit expression for the renormalized stress-energy tensor 
T a/ 3 in 4D (even for spherically-symmetric spacetimes). This makes it hard to construct 
the semiclassical BH geometry and evaluate its surface gravity. Fortunately the situation is 
much simpler in the two-dimensional (2D) framework, wherein is explicitly known for 
a generic background metric. This motivates us to address this issue — the mass- dependent 
semiclassical correction to the Hawking temperature and outflux — in the 2D framework. 

About two decades ago Callan, Giddings, Harvey and Strominger (CGHS) [2] introduced 
a formalism of 2D dilaton gravity in which the metric is coupled to a dilaton field and 
to a large number N of identical massless scalar fields. In this 2D framework T a p is known 
explicitly, allowing one to translate semiclassical dynamics into a closed system of partial 
differential equations (PDEs) [2]. Although the exact solution to these PDEs is not known 
explicitly, it is possible to explore these solutions numerically, and also through certain 
analytical approximations (see below) [3]. 

The purpose of this paper is to explore static BH solutions [I] of the 2D semiclassical 
CGHS model, 1 and to find out how the surface gravity (and hence also temperature and 
outflux) changes with the BH mass M, due to the semiclassical backreation on the metric, in 
the domain of large mass (M ^> m p ). We shall address this problem here both numerically 
and analytically. Note that in the static case the CGHS field equations reduce to ordi- 
nary differential equations (ODEs), which drastically simplifies their numerical (as well as 
approximate-analytical) solution. At the analytical side, since for large M the semiclassical 
BH geometry is well-approximated (locally) by the corresponding classical solution, we shall 
consider the semiclassical solution as a small deviation from the classical one, and treat this 
deviation by linear perturbation analysis. As it turns out, the overall magnitude of this per- 
turbation scales as 1/M. Not surprisingly, the perturbation analysis yields a semiclassical 

1 From the physical view-point, a static semiclassical BH solution should be viewd as the (somewhat hypo- 
thetical) situation in which a constant influx arrives from past null infinity, and exactly compensates for 
the Hawking outflux. 
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correction to k (as well as to T# and the outflux) which scales as 1/M too. We calculate 
this correction analytically, and then confirm it numerically. 

Note that there is a remarkable difference between 2D and 4D BHs, already at the classical 
level: Whereas in 4D the surface gravity k scales as 1/M, in 2D it is actually independent 
of the BH mass. As a consequence, in a semiclassical macroscopic BH the (leading-order) 
Hawking temperature Th is oc 1/M in 4D, but constant (i.e. independent of M) in the 
2D framework. Correspondingly the outflux (the outgoing component of T Qj g) at infinity is 
oc M~ 2 in 4D but constant in 2D. The semiclassical correction to the background geometry 
of the 2D BH modifies these constant values of k, Tjj and outflux, and the modification 
in all three quantities scales as 1/M at the leading order. It is this oc 1/M leading-order 
semiclassical effect which we explore in this paper — both theoretically (through perturbation 
analysis) and numerically. 

We point out that a similar phenomenon also occurs in the case of a dynamical, evap- 
orating, 2D CGHS BH. At the large-mass limit, a 2D BH evaporates at a constant rate 
M (owing to the M- independence of k in the classical CGHS solution). However, due to 
the backreaction of the semiclassical T a ^ on the geometry, there is a finite-mass correction 
to the Hawking outflux (and hence to M), which again scales as 1/M. 2 This correction 
term for a 2D evaporating BH was calculated analytically [5], and also confirmed numeri- 
cally [HI [7]. It is remarkable that the leading-order (oc 1/M) finite- mass correction to the 
(otherwise-constant) outgoing component of T a p is found to be exactly the same in the static 
and evaporating cases. We further comment on this observation in the Discussion section. 

In the next section we briefly review the CGHS semiclassical model (as well as its classical 
limit). We restrict attention to static solutions which are regular at the horizon, and explore 
the asymptotic behavior of the semiclassical geometry at the horizon and at infinity. We 
then turn in section 3 to address the Hawking temperature of such static BH's, as well as the 
Hawking outflux (and influx) at infinity, taking into account the deviation of the background 
geometry from its classical counterpart. Then in Sec. 4 we analyze the large-mass leading- 
order (oc 1/M) semiclassical correction to the geometry, treating it as a linear perturbation. 
This in turn yields the oc 1/M correction to k, and hence to the BH temperature T# and 

2 In the evaporating case, this "M" actually refers to the Bondi mass, namely the remaining BH mass as 
seen by an observer at FNI. 



4 



the Hawking outflux. These analytical results are verified numerically in Sec. 5. Finally in 
Sec. 6 we briefly discuss our results, and compare them to the corresponding case of 2D 
evaporating BH. 



II. THE CGHS MODEL 

The CGHS model [2] consists of a two-dimensional metric g a p coupled to a dilaton </> and 
to a large number iV ^> 1 of identical minimally-coupled, massless, scalar fields fi. It is 
convenient to express the metric in the conformal form, 

9u v 2 ^ ^ ' 9uu 9w 0. (2-1) 
At the classical level, the action then takes the form 

N 



(2.2) 



1 f 1 

Sdassic = - dudv e~ 20 (-2p ;UV + 4<f) ;U (f) >v - X 2 e 2p ) — - /<,„/»,» 

71 L t=i 

where A is a constant. At the semiclassical level, the trace anomaly contributes an effective 
term [2] 

N f 

Strace = / P ,uP ,vdud,V , 

leading to the overall semiclassical action 

Ssc S c l ass i c Strace- (^"^) 

Variation of S sc with respect to the scalar fields /j yields the standard wave equation 
fi,uv = 0, and throughout the paper we shall consider the trivial solution /, = 0. Variation 
of (ft yields a certain nonlinear hyperbolic equation, and the variation of the metric yields 
one hyperbolic equation and two additional constraint equations. Overall, there are two 
evolution equations and two constraint equations (which will be presented shortly). 



The coordinate transformations which preserve the double- null structure (2.1) of the line 
element are of the form u — > u(u), v — > v(v). They transform p into 

1 du 1 dv 
p = p In- l n ^~- 

2 du 2 dv 

Before analyzing the semiclassical dynamics, it will be useful to consider the (much sim- 
pler) classical system, construct its general static solution, and explore its asymptotic prop- 
erties. The insights gained from the classical system will in turn facilitate our analysis of 
the more complex semiclassical dynamics. 



A. Classical equations 



The classical field equations are obtained by varying S c i assic with respect to the dilaton 
and the metric. These equations are much simplified by the fields redefinition 

R = e~ 2 *, S = 2{p-<j>). 

The evolution equations then take the form 

R :UV = -e s , S :UV = 0. (2.4) 

In addition there are two constraint equations, which (after substituting f\ = 0) read 

R,uS iM ~ R,uu = , R iV S >v — R iVV = 0. (2-5) 

Note that 

p= (S-\nR)/2. 

In a coordinate transformation u — > u(u),v — > v(v), R is unchanged but S transforms 
according to 

~ „ , du , dv 

5 = 5 -In— -In— . 2.6 

du dv 

The form of the evolution and constraint equations is unchanged by such a transformation. 

The general solution of this system of classical field equations (evolution+constraints) is 
known to be uniquely described by a one-parameter family of solutions (up to coordinate 
transformation), which in the so-called Eddington-like coordinates take the form 

R = M + e v ~ u , S = v-u, (2.7) 

where M is an arbitrary constant. For M > (which we shall assume throughout) these 
solutions describe a black hole (BH) with ADM mass M. 

Throughout this paper we shall restrict attention to static (classical and semiclassical) 
solutions, namely, solutions which only depend on the spatial variable x = v — u. All field 
equations then reduce to ordinary differential equations (ODEs), which drastically simplifies 
the analysis. The classical evolution equations then read 

R" = e s , S" = 0, (2.8) 
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and the two constraint equations reduce to a single one: 



R'S' - R" = 0, (2.9) 



where a prime denotes differentiation with respect to x. The general classical solution (2.7) 
is manifestly static, and we re-write it as 

R = e x + M, S = x. (2.10) 

In this classical vacuum solution 

p = -[x-\n(e x + M)}. 

Note that p vanishes at the limit x — > oo (corresponding to spacelike as well as null infinity), 
implying that the classical BH spacetime is asymptotically flat. 

The other asymptotic limit, x — > — oo, corresponds to approaching the black- hole horizon. 
At that limit R decays exponentially to M, whereas S — > — oo. 3 

To facilitate the semiclassical analysis below, it will be useful to explore the general 



(static) solution of the classical evolution system (2.8), while relaxing the constraint equation 



( 2.9[ ). We shall refer to this more general solution of Eqs. (2.8) as the (static) classical 



flux-carrying solution — to be distinguished from the strict vacuum solution (2.10) (which 
satisfies the constraint equations as well). 4 As we shall see below, in the two important 
asymptotic regions — horizon and infinity — the semiclassical solutions are well approximated 
by certain flux-carrying classical solutions (which mimic the semiclassical fluxes in these two 
asymptotic regions). 



The integration of the evolution equations (2.8) is straightforward, and one obtains the 
general solution 

R = a - 2 e ax+b + cx + m, S = ax + b } (2.11) 

with four arbitrary constants a,b,c,fh (as expected for two second-order ODEs). However, 
by a simple, linear coordinate transformation of the form u — > u'(u), v — > v'(v) one can get 



3 The divergence of x and S (and p) at the horizon merely represents a coordinate singularity. To remove this 
singularity one may transform from the Eddington-like coordinates u, v to new, Kruskal-like coordinates 
U = -e~ u and V = e v . 

4 More generally a classical (non-static) flux-carrying solution will be defined to be a solution of the 



evolution equations (2.4 1 which does not necessarily satisfy the constraint equations (2.5) 
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rid of two parameters, say a and b, 5 bringing the solution to its canonical form: 



R = e x + cx + m, S 



x. 



(2.12) 



Thus, the static classical flux-carrying solution is actually a two-parameter family, in which 
c represents the flux magnitude (as we shortly discuss), and the other parameter m is 



reminiscent of the mass parameter [compare to Eq. (2.10)]. 

For any configuration R(x), S(x) we define the flux T(x) to be the differential expression 



at the left-hand side of the constraint equation (2.9), namely T = R'S' — R". In a flux- 



carrying solution, the evolution equation (2.8) guarantees that T = const, as one can easily 



verify . Expressing the flux-carrying solution in its canonical form (2.12), one readily finds 
that T = c. 



The spacetime's geometry is fully dictated by p(x), which may be expressed as 

1 



P 



In [l + (cx + m)e x ] . 



(2.13) 



At spacelike infinity (x — > oo), p oc (cx + m)e x . In particular this implies that 



(x ->■ + 



oo, 



(2.14) 



with all three quantities at the left-hand side decaying exponentially in x. Thus, as long 
as the strict spacetime metric is concerned, the classical solutions (vaccum as well as flux- 
carrying) are all asymptotically-flat: the metric approaches 2D Minkowski (p = 0), and the 
curvature asymptotically vanishes. 

Next we consider the asymptotic behavior of the flux-carrying solutions at the horizon 
limit (x — > — oo). It will be instructive to deduce this asymptotic behavior directly from the 



evolution equations (2.8) (rather than by just substituting x — > — oo in the explicit solution), 



5 The coordinate transformations which preserve the double-null form of the metric, as well as staticity 
(namly, dependence on x = v — u only) must be of the linear form u = au + /3 U , v = av + fi v (such that 



x 



ax + (3, where j3 = j3 v — j3 u ). By appropriate choice of a and /? one can get rid of a and b in Eq. 



(2.11 1, which also modifies c and rh. [Note that in this coordinate transformation R is unchanged and S 



changes into S = S — 2 In a: However, to address the transition from (2.111 to (2.12| one needs to consider 
the functional form of R(x) and 5*(i).] 
6 More generally, for any field configuration R(u,v), S(u,v) one may define the fluxes T uu and T vv to be 



the corresponding two differential expressions in the constraint equations (2.5). If the evolution equation 



(2.4) is imposed, it implies that d v T uu — d u T vv = 0. Then, if the solution is further assumed to be static, 



T U u = T vv = T and the latter must be a constant.) 
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because the same procedure will later serve us in the semiclassical case, wherein we shall 
analyze the near-horizon asymptotic behavior of the semiclassical evolution equations (see 
below) . 

The equation 5"' = dictates S' = const, yielding S — > — oo at the horizon. This, in 
turn, implies that the source term e s at the right-hand side of the evolution equation for 
R dies out, leading to R" = and hence to a linear asymptotic behavior of R(x) at the 
horizon. This is of course confirmed by the linear asymptotic form of the explicit solution 



(2.12) as x — >■ — oo. 

The Ricci scalar R® is given by —4e~ 2p p". At the limit x — > — oo one finds (for c ^ 0) 

R a a S -— e~ x -> oo. (2.15) 

[We also note that the scalar (V0) 2 too diverges like oc (c/x)e~ x as x — > — oo.] Thus, all 
flux-carrying solutions with c ^ develop a curvature singularity at the "horizon" limit 
x — > — oo. The only regular (static) BH solutions are those with c = (and M > 0), namely 
the true vacuum BH solutions. The c = solutions are perfectly regular at the horizon, and 
can be extended to the interior domain R < M in the usual procedure, by transformation 
to Kruskal-like coordinates (see footnote [3]) . 

For later convenience we summarize here the (classical) x — > — oo asymptotic behavior 
as follows: There is a one-parameter regular class which behaves as 

R = m } R' = 0(e x ) , S = x; (2.16) 

And there also exists a larger, two-parameter class, which behaves as 

R^cx + const , R' S c , S = x, (2.17) 

which however admits a curvature singularity rather than a regular horizon. 

B. Semiclassical equations 

At the semiclassical level the field equations get quantum correction terms coming from 
Strace- The semiclassical evolution equations become 

R,uv = ~e s - Kp iUV , S tUV = Kp iUV /R, (2.18) 
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where K = N/12. The two constraint equations now take the form 



R,ww R,w^,w K \P t wiv P yj z w (w)J — T ww , (2.19) 

where hereafter w will stand for either u or v, and by T ww we refer to the ww component 
of the renormalized stress-energy tensor. The functions z u (u),z v (v) carry the information 
about the quantum state, and should thus be determined by initial or boundary conditions. 

In the static model, where both R and S depend on x = v—u only, the evolution equations 
reduce to the ODEs 

R" = e s -Kp\ S" = Kp"/R. (2.20) 
The constraint equations again reduce to a single ODE: 

R'S' -R" = K [p" - p' 2 + z]=f, (2.21) 



where now T = T vv = T uu . Note that in the static case the w = u version of Eq. (2.19) 



implies z = z u (u), and similarly its w — v version implies z = z v (v), hence in Eq. (2.21) z 



must be a constant (it will be determined below by the regularity condition at the horizon). 
1. Far-field asymptotic behavior 

The asymptotic behavior of the semiclassical solution at infinity turns out to be fairly 
similar to the (flux-carrying) classical solution. Although we cannot prove this rigorously, 
this conclusion is easily derived from a simple iterative procedure, starting from the classical 
solution. Since in the latter p" decays exponentially in x, we can drop p" in the evolution 



system (2.20). The latter then reduces to the classical evolution system (2.8), which we 
have already integrated in the previous subsection. We find that the far-field semiclassical 
asymptotic behavior again takes the general form 

R^ e x + cx + m, S^x (x -> +oo) (2.22) 
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(with exponentially-small corrections), where c and m are free parameters 7 . In particular, 
in the semiclassical case p(x) again satisfies the asymptotic behavior 



p,p',p"^0, 



[x — > +00) 



(2.23) 



with all three quantities decaying exponentially (which in retrospect guarantees the consis- 
tency of this far-field approximation). 



Note that when the asymptotic expression (2.22) is substituted in the semiclassical con- 
straint equation (2.21), it yields 8 

z = c/K. 



2. Near-horizon asymptotic behavior 



We turn now to analyze the semiclassical asymptotic behavior at the horizon limit x — > 
—00. To this end, it is useful to represent the evolution equations in their "canonical" form 
(where R" and S" are explicitly expressed in terms of R and S). Substituting p = (S— In R)/2 



in Eq. (2.18) and then extracting R" and S" , one obtains 



R" 



and 



S" 



S 2R-K_ _ 
2(R-K) 

K 



K 



+ R 



2R (R - K) 
K 



/2 



(2.24) 



(2.25) 



2R{R-K) ' " 2R 2 (R-K)' 
To analyze these equations we again apply an iterative procedure, starting from the (regular) 
classical near-horizon asymptotic behavior (2.16). Since both R! and e s decay exponentially 



on approaching the horizon (x — > —00), in the two evolution equations ( 2.24|2.25 ) the right- 
hand side decays like e x too, namely R" = S" = 0(e x ). Integrating these equations, one 
finds that both R' and S' approach finite values at the horizon: Rl — y const = Ri,S' — > 
const = S\, where Ri and Si are yet arbitrary. However, one again finds that as long as 
Ri 7^ 0, the Ricci scalar diverges like oc Rie~ x — > 00 at the horizon [in full analogy with 



7 Just like in the classical case, the general solution of the evolution equations admits four free parameters, 
but by applying a gauge transformation we get rid of two parameters (see footnote [5] above) and obtain 



the canonical form (2.22) 



By iteration one finds that the large-x corrections to Eq. (2.22) are 0(xe x ) for R and 0(xe 2x ) for S 



(with a pre-factor oc Kc in both cases), which at x — > 00 yield no contribution to Eq. (2.21 ) 
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Eq. (2.15) above]. For obvious physical reasons we restrict our attention here to horizon- 
regular semiclassical BH solutions, and therefore we set R\ = 0. Doing so, we obtain the 
general regular near-horizon asymptotic behavior: 

R ^ const = m , S = S + S x x, (2.26) 

with 0(e x ) deviations. The derivatives of R and S satisfy 

S' -»■ const = St (2.27) 

and 

R> = R" = S" = 0(e x ). (2.28) 
Setting p' = (S' — R'/R)/2, one finds at the horizon limit: 

S 

p ->■ const = y , p" ->■ 0. (ar ->■ -oo) (2.29) 

Recall also that in the classical solution S\ = 1. 

A few remarks are in order here: (i) In principle we could easily bring the value of Si 
to unity by a coordinate transformation, as we have done in the classical solution. Such 
a transformation, however, would multiply S'(x) by a global parameter (the constant 1/at 
in the terminology of footnote |5l) . But we have already fixed this constant by demanding 



S' = 1 at spacelike infinity, as demonstrated in Eq. (2.22) (see also footnote [7]) ; hence in the 
semiclassical case we no longer have the freedom to gauge out the near-horizon parameter S% . 
(ii) m is a natural extension of the classical mass parameter m to the semiclassical case, (iii) 
The above near-horizon analysis is restricted to the case m > K: Then, when R approaches 



m at the horizon, the terms R — K in the denominators in Eqs. ( 2.24|2.25 ) remain bounded 
away from zero, (iv) In fact, throughout the rest of this paper we are primarily interested 
in the macroscopic case, namely rh ^> K. In this case, the solution is well approximated 



by the classical solution. In the latter S\ = S' — 1, as seen in e.g. Eq. (2.12). For this 
reason, in a macroscopic BH one expects that S\ ~ 1 (this will be shown more explicitly 
in the next section, in which the deviation of Si from unity is explicitly calculated in the 
large- mass case and shown to vanish at the macroscopic limit). In particular, it follow that 
Si is positive (at least in the large-mass case). 
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We also point out that the asymptotic behavior ( 2.27|2.28 ) guarantees that R' and e s 
indeed decay exponentially on approaching the horizon, which in retrospect justifies our near- 
horizon approximation [in which we have regarded the right-hand sides of Eqs. ( 2.24|2.25 ) 
as negligible]. 



3. Concluding remark 



In the previous section we asserted that the flux-carrying classical solution well approx- 
imates the semiclassical solution at both asymptotic boundaries, namely at the far-field 
region and near the horizon. We are now in a position to verify and further clarify this 



statement, by comparing the exact classical solution (2.12) to the semiclassical asymptotic 
solutions at x — > oo and x — > — oo. 



In the far-field region, Eq. (2.22) conforms with (2.12), with the mapping c — )■ c,m — >■ m 



(and, as will become clear in the next section, c ^ in the semiclassical case — which is the 



essence of the Hawking effect). Near the horizon, Eq. (2.26) conforms with Eq. (2.12), but 
this time with the mapping c — > 0, m — > fh. 

So here is a key difference between the classical and semiclassical solutions: In the classical 



case, the exact flux-carrying solution (2.12) holds all the way from infinity to the horizon, 



with a single set of parameters c, m. Thus, regularity at the horizon (c = 0) automatically 
guarantees vanishing flux at infinity (and everywhere). On the contrary, in the semiclassical 



case, although Eq. (2.12) is a good approximation both at the near-horizon and far-field 



regions, it fails to approximate the solution in between, and the consequence is that the 
free parameters c, m take different values at the two asymptotic regions. In particular, 
the difference in c (namely c = at the horizon but c ^ in the far-field region) is a 
manifestation of the Hawking effect, as we further discuss in the next section. 

Notational remark: Throughout the rest of this paper, we shall replace the symbol fh by 
m for notational simplicity. Thus, by the mass m of a semiclassical BH, we shall specifically 
refer to the value of R at the horizon. 
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III. HAWKING TEMPERATURE AND THE SEMICLASSICAL FLUXES 



A. Surface gravity and Hawking temperature 

The Hawking temperature T# of a static BH is directly dictated by the horizon's surface- 
gravity k. When the 2D metric is expressed in the Schwarzschild-like form 

ds 2 = -F(r)dt 2 + [1/F(r)]dr 2 , 

the surface gravity is given by the standard expression 

1 fdF\ , . 

K= 2\^) ■ (3 - 1} 

^ \ / (horizon) 

In the double-null metric we use here, k takes the simple form 

K = 2f ''(horizon)- (3-2) 

[To verify this, one may start from the above Schwarzschild-like line element, define r*(r) 
through dr/dr* = F(r), and then introduce the null coordinates v = t + r*, u = t — r*, 



recovering the metric (2.1) with e 2p = F(r). Noting that x = 2r*, one obtains Eq. (3.2) 



In our 2D semiclassical model, setting 2p' = S' — E! /R and recalling the near-horizon 



asymptotic behavior ( 2.27[2.28 ), we obtain 

k = S v (3.3) 



Note that p' in Eq. (3.2) is in fact gauge-dependent, and the right-hand side of Eq. (3.2) 



is to be evaluated using the "asymptotically-Minkowskian gauge" , namely the specific gauge 



in which p vanishes at infinity [which is indeed the gauge we use here, cf. Eq. (2.23)]. For 
later convenience we also express k in a gauge-invariant manner: 



2P(-oo) _ ^(-oo) 

oi q> 



(3.4) 



where ^(±00) = S'(x —> ±00) etc.: In a gauge transformation x = ax + (3 (cf. footnote 
[5j), the difference ds (or dp) associated with two adjacent spacetime points is unchanged, 
whereas dx changes by the factor a, which however cancels out between the numerator and 
denominator. 
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Once k is obtained, the Hawking temperature Th is given by [8] 

T H = k/2tv. (3.5) 

This relation may be derived straightforwardly by analytically extending the 2D metric (e.g. 
in the above Schwarzschild-like form) into the euclidean domain, just like Hawking's original 
derivation pQ in the 4D case. 

At the classical limit k = S± = 1. It is remarkable that unlike the 4D Schwarzschild BH, 
wherein k oc 1/m, in the classical 2D case n is independent of m — and so is the Hawking 
temperature. 



B. Semiclassical fluxes at infinity 



As was discussed above, in a static solution the semiclassical outflux and influx are equal 
at each point: T vv = T uu = T(x). One of our main goals in this paper is to calculate the 
value of T at infinity, which we shall denote T^. 

Equations (2.21 ) and (2.23) yield = Kz. The parameter z, in turn, may be determined 



from the asymptotic behavior at the horizon: Recalling Eqs. ( 2.27[2.28 ), the left-hand side 



of Eq. (2.21) vanishes, which by virtue of Eqs. (2.29) and (3.3) implies 

si 



The fluxes at infinity are therefore 



T 



K 



K 

T 



K- 



(3.6) 



4 4 

Alternatively the outflux at asymptotic infinity may be obtained from the Hawking 



temperature Th- Owing to asymptotic flatness of the 2D metric [cf. Eq. (2.23)], we may 
apply at asymptotic infinity the standard rules of flat-space 2D physics 9 : For each of the 
iV scalar fields /j, the outflux is related to the temperature through the 2D analog of the 



Stefan-Boltzmann law: = (ir 2 /12)T H . Multiplying by iV = 12K and recalling Eq. (3.5), 



one recovers Eq. (3.6) 



9 Note that in the CGHS 2D model the scalar fields /, do not admit any scattering due to curvature, hence 
there is no "gray factor" in the Hawking outflux (unlike the 4D case). 
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Thus, all that is needed for the determination of is Si, namely the horizon value of S'. 
At the macroscopic limit (m — >■ oo), one can employ the classical solution for this purpose, 
yielding Si = 1 and hence = K/4. For a finite m there will be a deviation from this 
classical value of Si, which will grow with decreasing mass. In the next section we shall 
calculate this deviation for large m, at the first order in m/K. 



IV. FINITE-MASS SEMICLASSICAL CORRECTIONS 



The spacetime outside of a macroscopic semiclassical BH is (in a local sense) well approx- 



imated by the classical solution. This is most easily seen in the evolution equation (2.25) for 



S, wherein the semiclassical correction is proportional to K. When substituting the classical 
expressions for R and S in the right-hand side, and taking the leading order in K/R, one 
recognizes that each of the two terms in the right-hand side is bounded by K/2R — which in 
turn is bounded, outside the BH, by ~ K/m. We may therefore expect that the semiclassical 



correction to S should be small whenever m 3> K. Inspecting the evolution equation (2.24) 



for R leads to a similar conclusion: The semiclassical contribution to R" is smaller than the 
classical contribution (e s ) by a small factor of order K/R < K/m. This conclusion — namely 
the smallness of the semiclassical correction to both R and S in the case m 3> K — was also 
verified in our numerical simulations (described in the next section). 

Thus, to investigate the small semiclassical correction in the case of a large-mass BH, it 
is useful to express the solution in the form 

R(x) = R c i(x) + 8R(x) , S(x) = S c [(x) + 8S(x), 



where R c i and S c i denote the classical solution (2.10). The semiclassical corrections 5R 



and 5S can then be treated as small perturbations. Linearizing the evolution equations 
( p34][2T25| in 5R and 5S (and in K), one obtains the inhomogeneous linear system 



5R" = 5Se s < 1 + 



K 



2R. 



d 



(R, 



R 



d 



and 



SS" 



K 



[(R'd? ~ Rcie 3 -} 



Substituing the explicit expressions (2.10) for R d and S c /, the system takes the form 

me x 



5R" = 5Se x + K- 



2(e x +mf 



(4.1) 



(4.2) 



(4.3) 
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and 



6S" 



Km 



(4.4) 



2 (e x + m) 

The integration of this system is straightforward, and the full explicit solution is presented 
in the Appendix. Here we shall only need the equation for SS, because the parameter k — Si 
(which determines the Hawking temperature and outffux) is entirely determined by 5S'(x). 



It is convenient that Eq. (4.2) for 5S" is decoupled from 5R. Its first integral yields 

ss> Km 



\2 ' 



(4.5) 



4 (e x + m) 

Notice that we have chosen the integration constant such that 5S' vanishes at x — > oo, to 



comply with the canonical form (2.22) of the asymptotic behavior at spacelike infinity. 



At the horizon limit we obtain 



S 1 = 1 + 8S'(x ->• -oo) = 1 + AT/ (4m). 



Thus, up to first order in K/m, 



k = 1 + K/(4m) 



(4.6) 



Substituting this in Eqs. (3.5) and (3.6) , we finally obtain 



T, 



ii 



for the Hawking temperature, and 



T 



2tt 



K 
T 



4m V m 



K „(K 
2m \m 



(4.7) 



(4- 



for the outflux. 

It is interesting to note that a similar result was previously obtained for the Hawking 
outflux emitted from a non-static, evaporating, CGHS semiclassical BH. In that case too, 



the outflux was found to take precisely the form (4.8), this time with m denoting the Bondi 



mass (namely the remaining BH mass, as measured at FNI as a function of u). That result, 
too, was obtained analytically [5] and then confirmed numerically P,[7j. Here we were able 
to get this first-order result with a much simpler calculation (both analytical and numerical), 
for a static BH model. 
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V. NUMERICAL VERIFICATION 



We carried out a numerical calculation to verify our analytical result (4.6) for the first 



order finite-mass correction to k. To this end we numerically integrated the static semiclas- 
sical evolution equations, namely the nonlinear ODE system ( 2.24|2.25 ). 

The semiclassical CGHS model admits a unique scaling law: [31 E] The field equations 
remain invariant under a multiplication of both K and R by the same arbitrary constant 
(while keeping p unchaged). It follows from this scaling law that the dependence of k on 
m and K (in the exact semiclassical solution) is only through the combination m/K. In 
our numerical simulations we calculated k as a function on m/K throughout the range 
100 < m/K < 500. 



Figure 1 



1.0025 



1.002 



---Numerical result 

First-order analytical result 

Zeroth-order expression 




1.0015 



K 



1.001 



1.0005 



1 

100 200 300 400 500 

m/K 



FIG. 1: The surface gravity k as a function of the BH "normalized mass" m/K. The dashed- 
dotted black curve displays our numerical results, which were calculated for 250 values of m/K 
throughout the range 100 < m/K < 500. The first-order theoretical prediction k = 1 + A;/ (4m) is 
the solid red curve. These two curves are actually visually indistingwishable. For comparison we 
also display the corresponding zeroth-order expression k = 1, the solid green line. 
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In the analytical derivation above we have imposed the regularity condition •R(_ DO ) — 
at the horizon. In addition we have worked with the gauge in which S^) = x, S'^ = 1, 
in accord with Eq. (2.22). In a numerical implementation this would mean that we would 
have to impose boundary conditions at both boundaries of the simulation, x — > — oo and 
x — > oo. This is somewhat inconvenient, though it can still be done, e.g. by iterations. 

However, to simplify the numerical procedure we have chosen a different approach: We 
can use a more convenient gauge for the simulation, and then calculate k using the gauge- 
invariant expression (3.4). We have chosen the gauge for the simulation to be S^-oo) = 
x, S^-oo) = 1 at the horizon, and integrated the ODE system ( 2.24|2.25 ) from the horizon 
towards infinity. The boundary conditions for R are taken to be R't-oo) — as well as 
R(-oo) = fa. In this way all our boundary conditions are specified at x — — oo. In this gauge 



it easily follows from Eq. (3.4) that k = 

The results of our numerical simulations are presented in Fig. 1. The graph demonstrates 
an excellent agreement between the numerical solution and our theoretical prediction, based 
on the first-order analysis. This figure also indicates very clearly that the first-order correc- 
tion is really needed: The numerical results deviate substatially from the zero-order outflux 
(the horizontal line k — 1, which is the infinite-mass limit). 

Finally we note that our numerical simulations also reveal a second-order correction term 
in k. This additional term is consistent (within the numerical accuracy) with 

3 (K x 1 
32 \m 

This of course implies a corresponding second-order correction in the Hawking temperature 
and the outflux. 



VI. DISCUSSION 

We have found analytically, and confirmed numerically, that in a large- mass (m ^> K), 
static, CGHS semiclassical BH the surface gravity k admits a small semiclassical correction 



oc K/m, presented in Eq. (4.6). Correspondingly, the Hawking temperature and outflux 



also get finite-mass corrections, as shown in Eqs. ( 4.7[4.8 ) 



As was already mentioned above, previous analytical [5] as well as numerical [SI [7] analyses 
of the Hawking radiation emitted from an evaporating CGHS BH revealed the presence of a 
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semiclassical correction term oc K/m in the outflux. Our analysis shows that this correction 
term is precisely the same (at order Kj m) in the static and evaporating cases. Namely, Eq. 



(4.8) holds in the evaporating well (provided that one interprets m as the Bondi 

mass) . 

In the evaporating case, one may be tempted to interpret the finite-mass correction 
to the outflux as an indication for the non-thermal character of the Hawking radiation. 
Such deviations from thermality might be important for certain aspects of the information 
puzzle Our analysis of the static case suggests a different interpretation of this oc K/m 
correction to the outflux: The backreaction of the renormalized stress-energy tensor T a p 
obviously modifies the background geometry, and as a consequence the BH's surface gravity 
k is modified too. This inevitably leads to a change in the Hawking temperature, through 
the standard relation Th = n/2ir (which should exactly hold, for a static BH, even in 
the semiclassical case). In turn, this correction in Tjj naturally leads to a corresponding 
correction in the outflux, through the standard quadratic relation between temperature and 
outflux (the 2D analog of the Stefan-Boltzmann law, which states that the outflux should 
be proportional to the square of the temperature). Indeed, this exact quadratic relation 
between outflux and Hawking temperature is guaranteed to hold in the 2D semiclassical 



BHs of the CGHS model, as demonstrated in Eqs. (3.5) and (3.6). This relation is also 



reflected in our explicit first-order results, Eqs. (4.7) and (4.8). Thus, in the static case the 
Hawking outflux remains precisely thermal (despite the semiclassical backreaction), even 
though the Hawking temperature gets a small semiclassical correction. 

Turning now to the evaporating case, since the first-order correction to the outflux pre- 



cisely agrees with the static result (4.8), it appears that this 0{K/m) correction merely 
reflects the change in the Hawking temperature due to backreaction (rather than devia- 
tions from thermality). Indeed one may anticipate small deviations from thermality in the 
dynamical, evaporating case; however, such deviations seemingly appear only at second or 
higher orders in K/m. 
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Appendix A: Linear semiclassical correction 



We display here the full solution to the linearized semiclassical equations (|4.3|4.4|): 

K 



5R 



5S 



In 



K 
4m 



e x + m 



m 



rn 



H (x — \n(e x + rn)) 

m 



+ x — In (e x + m) 



e x + m 

The four integration constants were chosen here such that 5S and SS' vanish at infinity (in 
accord with our gauge condition), and 5R, 5R' vanish at the horizon. [By differentiating the 
last equation one recovers Eq. ( 4.5[ ) for 5S'.) 
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